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THE  CALCULATIONS  OF  TWO  COMPRESSIBILITY 
MECHANISMS  IN  EQUATIONS  OF  SOIL  DYNAMICS 


G.  M.  Lyakhov  and  N.  I.  Polyakova 


(Moscow) 


A  system  of  equations  is  proposed  which 
determines  the  unidirectional  movement  of  the 
soil  with  the  simultaneous  calculation  of  its 
compressibility  as  a  three-component  medium 
and  also  as  the  compressibility  of  the  soil 
skeleton.  The  soil  characteristics  and  the 
values  of  the  pressures  at  which  it  is  neces¬ 
sary  to  calculate  both  compressibility  mechan¬ 
isms  were  determined  experimentally. 


1.  Soil  is  an  aggregate  of  solid  mineral  particles  (grains) 
which  form  its  skeleton,  and  pore  spaces  between  these  particles 
which  are  filled  with  water  and  air.  During  the  action  of  a 
dynamic  load  deformation  of  the  skeleton  occurs  along  with  the 
simultaneous  compression  of  the  water,  air  and  solid  particles. 
In  accordance  with  this  two  mechanisms  of  compressibility  and 
two  mechanisms  of  stresses  need  to  be  examined  [1,  33* 


The  first  mechanism  corresponds  to  the  compressibility  of  the 
soil  skeleton.  The  deformation  is  associated  with  the  destruction 
of  the  bonds  between  the  solid  particles  and  their  displacement 
relative  to  one  another.  Experiments  indicate  that  there  is  stress 
in  soils  as  well  as  in  solids,  and  depending  on  the  orientation 
the  state  of  stress  is  determined  by  the  tensor  value. 
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The  second  mechanism  corresponds  to  the  compressibility  of  ti. 
soil  as  a  three-component  medium,  i.e.,  a  mixture  of  solid  particles 
water  and  air.  The  deformation  of  the  soil  is  connected  with  the 
volume  deformation  of  gaseous,  liquid  and  solid  particles.  Prom 
experiments  it  follows  that  the  state  of  stress,  just  as  in  liquids, 
is  determined  by  the  scalar  gravity  —  the  pressure. 

During  the  dynamic  processes  one  of  the  stresses  frequently 
turns  out  to  be  substantially  smaller  and  it  can  be  ignored.  At 
a  stress  where  the  volume  deformation  of  the  soil  e  is  substantially 
less  than  the  allowance  for  the  volume  of  air  in  the  soil,  the 
compressibility  is  determined  generally  by  the  compressibility  of 
the  skeleton.  When  e  >  a^,  the  compressibility  of  the  soil  is 
determined  generally  by  its  compressibility  as  three-component 
medium.  When  e  approaches  a^,  it  is  necessary  to  take  into  account 
both  compressibility  mechanisms.  For  example,  in  water  saturated 
soil  with  a1  =  0.01  under  a  load  of  several  atmospheres  and  in  un¬ 
saturated  soil  with  a1  =  0.3  under  loads  of  several  tens  of  at¬ 
mospheres  . 

Let  us  find  the  equations  which  determine  the  soil  movement 
taking  into  account  both  compressibility  mechanisms.  Let  us  first 
examine  the  soil  properties  as  a  three-component  medium.  During 
short-term  dynamic  processes  the  value  of  the  components  in  each 
of  the  soil  particles  according  to  the  mass  remains  unchanged, 
but  by  volume  is  changed.  Let  us  designate  o^,  c*2,  a^,  as  the 
value  by  volume,  p^,  p2,  p^,  —  density,  and  c^,  c2,  c^,  —  the 
speed  of  sound,  respectively,  in  the  gaseous,  liquid  and  solid 
components  at  atmospheric  pressure. 

Then,  the  soil  density  pg  and  bulk  density  (volume  weight)  of 
the  skeleton  Yq  can  be  determined  by  the  expressions: 

Po  —  ai  Pi  +  «s  Pi  +  °a  P».  To  —  ®s  P».  (1*1) 
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In  this  case 


a,  -r  a,  -r  a,*  1. 


Based  on  the  assumption  that  with  applied  stress  p  the  pres¬ 
sure  or.  all  the  components  is  also  equal  to  p,  and  that  the  density 
of  each  component  car.  be  determined  by  law,  corresponding  to  its 
compression  in  the  free  state,  the  equation  of  the  state  of  the 
soil  [1]  as  a  three-component  medium  was  obtained: 


where  rt,,  * ... ,  * ..  —  indices  cf  tne  lse.ntropes  of  the  air,  water  and 

-L  <-  J 

quartz,  accepted  as  equal  to  1;  4;  3* 


The  unloading  of  the  soil  as  a  mixture  proceeds  according  to 
this  same  equation.  In  Fig.  1  the  graph  of  the  relationship  p  - 
=  p(e),  corresponding  to  (1,  2),  is  represented  by  curve  I.  De¬ 
formation  e  is  associated  with  the  density  of  the  soil  p  and  the 
specific  volume  V  by  the  expression 


-  YzlYs. 

P  V' 


Let  us  examine  the  properties  of  the  soil  as  a  solid.  It 
is  assumed  that  in  the  case  cf  or.e-dimens ional  planes  and  spherical 
movements,  the  soil  properties  are  determined  by  two  relationships 
known  from  the  experiment  -  the  lav;  of  cubic  compression  and  the 
condition  of  plasticity 


«=*«.  + 2  $  (p,  p0), 

/  0>i.  °3)  =  0. 


r  TD-  i  i'i  —  2  3—  7  ^  9  —  7  w 
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(l.*) 
( 1  • r;  ) 


Here  and  o0  -  the  stresses  acting  in  the  direction  of  the  ap¬ 
plication  of  the  external  load  as  well  as  in  a  perpendicular  di¬ 
rection,  <3  -  average  normal  stress,  p  -  current  stress,  and  p„.  - 

ill 

previously  attained  maximum  soil  density. 

In  the  first  approximation  with  stresses  not  too  great 

a>  =  #t  0  (°l)  «!•  (  1  .  C  ) 

In  many  soils  =  const.  From  here  (1.^4)  it  can  be  express¬ 
ed  in  the  form 

3|  (1  +  2  «,)  =  4  (p,  p») 

or 

«i  =  ?  («i.  •«).  ( i  ■  7 ) 

when  t  —  the  maximum  deformation  previously  proven  by  the  soil, 
m 

Let  us  examine  the  movement  of  the  medium  based  on  the  com¬ 
pressibility  mechanisms.  During  quick-acting  processes  the  water 
and  air  are  not  expelled  from  the  soil;  they  appear  to  be  forced 
into  the  pores.  Therefore,  it  can  be  assumed  that  each  particle 


U 


of  the  soil  making  up  all  three  components,  acts  as  an  indivlu:. 
unit  with  a  velocity  u  ■  u(x,  t)  coinciding  with  the  direction  of 
the  external  loan  (x  -  spatial  coordinate,  t  -  time).  The  movement 
of  the  medium  is  therefore  associated  with  the  overall  stress  in 
the  soil  P,  equal  to  the  sum  of  pressures  p  in  the  soil  as  a  mix¬ 
ture,  and  the  stresses  in  the  skeleton  o^,  acting  in  the  direction 
of  the  movement.  From  here  it  follows  that  the  unidimensional 
movement  is  determined  by  a  system  of  equations: 


Ajl 

dt 


du 


JU..O.  „.£  + 

x  dt 


d 

dx 


(p  —  ei)  — 


a.  —  a, 


L. 


0. 


(1.6) 


With  planar  and  spherical  symmetry,  v  is  equal  to  C  ar.d  2, 
respectively,  hquaticns  (l.d),  (1.2)  in  conjunction  with  (1.7) 
and  (1.5)  or  with  (l.e)  ana  (1.5))  form  a  closed  system  of  five 
equations  with  five  functions:  o^,  c^,  p,  p,  u.  This  system 
pertains  to  the  increment  cycle  of  the  external  lead,  i.c.,  to 
the  periou  of  increase  in  the  deformation. 

During  unloading  the  soil  density  as  a  mixture  is  lowered 
according  to  the  law  of  loading,  and  the  deformations  in  the 
skeleton  change  according  to  a  law  different  from  that  of  the  law 
of  leading.  During  unloading  =  o^e,  e )  ^ ^  ,  where  -  the 

maximum  deformation  attained  during  loading.  In  this  case 

e'"  >  *<n,  «i  (s,  '-)  >  0,  (  *,  (1.9) 

The  Dependences  a^(e)  during  loading  and  unloading  are  represented 
by  curves  2  and  3  in  Fig.  1.  Curves  4  and  5  pertain  to  the  sum 
P  =  p  -  o-p  during  loading  and  unloading,  respectively. 

Unloading  in  various  soils  can  proceed  differently.  In  soils 
where  the  rate  of  water  loss  is  small,  the  water  and  air  just  as 


during  loading  remains  In  the  restrained  state.  Each  soil  partial, 
expands  as  a  single  unit,  and  maintains  the  common  rate  of  all  the 
components.  The  movement  is  determined  by  this  same  system  of 
equations  also  by  the  fact  that  during  loading  only  the  stress  in 
the  skeleton  is  determined  by  the  equation  -  d^e,  e*'1'1). 

With  a  reduction  in  the  external  load  when  e  =  (Fig.  1) 

the  stress  in  the  skeleton  falls  to  zero.  All  of  the  load  ir.  this 

( ? ) 

case  is  taken  only  by  the  soil  as  a  mixture.  When  e  <  e  J  the 
expansion  of  the  soil  as  a  mixture  causes  the  formation  of  tensile 
stresses  in  the  skeleton.  In  this  case 


P  >  I  «!  I- 


The  system  of  equations  which  determines  the  movement  of  the 
soil,  in  this  case,  is  valid.  With  a  reduction  of  the  external 
load,  the  tensile  stresses  increase,  and  the  deformation  is  reduced. 
At  the  moment  the  external  load  drops  to  zero,  when  the  deformation 
of  the  soil  is  equal  to  £  (Fig.  1),  pressure  p  is  positive, 
and  the  stress  is  negative.  In  this  case 

P=  I  «i  I , 

e  ^  ^  corresponds  to  the  residual  deformation  of  the  soil.  The 

greater  the  content  the  water  in  the  pores,  the  less  intense  the 

(3) 

expansion  of  the  soil  e  . 


Subsequently,  after  the  termination  of  the  action  of  the 
dynamic  load,  the  water  and  air  gradually  begin  to  move  to  the 
side  of  the  free  surface,  and  the  pressure  in  the  soil  as  a  mixture 
falls  to  zero.  This  process  is  not  directly  connected  with  the 

propagation  of  the  wave.  In  soils  with  a  high  rate  of  water  loss 

( 2 ) 

during  unloading  when  e  <  e  ,  displacement  of  the  liquid  and 
gaseous  components  can  take  place  relative  to  the  skeleton.  In 
this  case  the  above  examined  system  of  equations  is  no  longer 
satisfied. 
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Certain  water  saturated  sands  pass  into  the  state  of  1 ’  'ue- 
f action  during  the  propagation  of  blast  waves.  This  associate 
with,  apparently,  the  arise  of  tensile  stresses  in  the  soil  after 
the  passage  of  the  shock  wave  front.  Similar  reasoning  was  offered 
by  V.  I.  Belokopytov. 

From  the  experiment  [3]  it  is  known  that  with  a  reduction  of 
the  content  of  air  ir.  the  water  saturated  sand,  the  probability 
cf  the  onset  of  liquefaction  is  increased.  The  following  explana¬ 
tion  can  be  given  for  this.  With  a  decrease  in  a,  curve  i  in  Fig.  1 
approaches  the  ordinate  axis.  Therefore,  during  similar  external 
leads  in  a  sand  with  a  smaller  a.,  the  values  of  p*'"^  will  be 

( i  \  1 

larger  and  o  smaller.  The  maximum  possible  value  cf  the  tensile 

1 1 ) 

forces  increases  with  an  increase  in  p'  .  In  this  way  with  a 
decrease  in  the  value  of  p  ^  ^  and  consequently  the  probability 
of  liquefaction  increase. 

After  liquefaction  sets  in,  the  movement  of  the  soil  is  deter¬ 
mined  by  the  equations  (1.2)  and  (1.8),  in  which  o-^  =  =  C.  The 

velocity  of  the  shock  wave  front  2  and  the  speed  cf  sound  c  in  the 
calculation  of  the  two  compressibility  mechanisms  are  determined 
by  the  equations: 


D  ,/  f  C-  , 

/^-8l)  (1.10) 

1  Po  ?-Po  V 

d  P 

where  p  is  associated  with  p  by  equation  (1.2). 

2.  Let  us  examine  the  results  of  the  experiments  for  the 
determination  cf  the  dynamic  compressibility  of  the  soil  and  the 
conditions  under  which  it  is  necessary  to  calculate  the  two  com¬ 
pressibility  mechanisms.  The  experiments  were  conducted  on  a 
sandy  soil  in  a  water  saturated  and  unsaturated  condition.  The 
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soil  was  placed  In  a  cylindrical  tank  with  cement  walls  and  be.. 

(the  diameter  of  the  tank  was  2  m,  the  height  1.85  m). 

The  volume  weight  (bulk  density)  of  the  skeleton  yn  = 

=  1*50-1.56  g/cm  .  The  amount  of  air  by  volume  in  the  unsaturated 
soil  *  0.25-0.37  (moisture  W  =  4-105),  and  in  the  water  saturated 

soil  =  0.02-0.04  (W  =  23-275)*  The  grain-zize  distribution  is 

given  below: 


Table 


\  MM  ! 

1 

:  1  |  1-0/, 

j  0,5—0,25 

0.25  0,1 

0.1-0.01 

<  0.01 

t  ...  1 

! 

1  is- is 

j  30—10 

|  15-25 

5-8 

2—4 

Here,  <$  —  diameter  of  the  particles  in  mm,  3  —  their  content  in 
percent . 

Waves  were  produced  by  a  blast  of  a  slab  charge  placed  on  the 
surface  of  the  soil  at  the  top  of  the  tank.  From  the  surface  the 
charge  disturbed  a  layer  of  soil  having  a  thickness 

A  =  *  C; 

2 

C  -  density  of  the  charge  (weight  in  kg  per  1  m  ), 

*  -  3  **/«£> 

The  parameters  of  the  waves  were  recorded  with  the  aid  of 
strain  gagos  whose  readings  were  recorded  on  loop  oscillographs. 

The  experiments  indicated  that  the  wave  is  a  shock  wave  when 

the  maximum  pressure  is  greater  than  p  .  In  water  saturated  and 

s  2 

uncatured  soils  p  is  equal  to  7-10  and  4-6  kg/cm  . 
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In  Fig.  2  graphs  of  the  dependence  of  the  wave  front  velcci 
D  and  the  maximum  pressure  Dm  on  the  relative  distance 


(2.1) 

C 

where  R  -  distance  from  the  site  of  the  (slab)  charge,  m;  C  - 

p 

density  of  the  slab  charge,  kg/m  ;  Curves  1,  2  pertain  to  the  water 

saturated  soil,  and  3-^  —  to  the  unsaturated  soil.  When  the  wave 

is  a  shock  wave,  D  and  coincide. 

m 


From  the  graphs  it  is  obvious  that  at  close  distances  to  the 
slab  charge,  the  velocities  in  the  water  saturated  soil  is  sub¬ 
stantially  higher  than  in  the  unsaturated  soil.  With  an  increase 
in  R°  the  differences  diminish. 

Figure  3  portrays  graphs  of  the  dependence  of  the  maximum 
pressure  on  R°  in  water  saturated  soil  (curve  1)  and  in  un¬ 
saturated  soils  (curve  2).  The  scattering  of  the  experimental 
points  is  great  which  is  associated  with  the  fluctuations  of 
a1  and  yQ  of  the  first  order  in  the  various  experiments.  However, 
all  the  points  pertaining  to  water  saturated  soil,  are  appreciably 
higher.  The  maximum  pressure  in  water  saturated  soil  is  reduced 
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significantly  more  slowly  with  distance.  The  found  dependence  of 
the  maximum  pressure  on  the  distance  can  be  presented  in  the  form 


(2.2) 


3. 


For  water  saturated  sand  «  =  55,  U  =  1.3,  and  for  unsaturated  sand 
h  =  12,  u  =  1.5. 

Let  us  construct  the  relationship  P-^  =  P^(e)  according  to  the 
experimental  data  of  P-^  and  D  in  the  range  of  pressures  where  the 
wave  is  a  shock  wave.  At  the  shock  wave  front  the  relationship  is 
satisfied 


(2.3) 


Here,  Pq  is  taken  to  be  equal  to  zero. 
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Graphs  =  P-^(e)  in  Fig.  ^  correspond  to  the  dynamic  com¬ 
pression  at  the  shock  wave  front  of  a  water  saturated  (curve  1) 
and  unsaturated  (curve  2)  soil.  They  were  plotted  in  accordance 
with  the  equation  (2.3)  based  on  the  experimental  curve  in 
Figs.  2  and  3-  Curves  3-3  express  the  relationship  p  =  p(e), 
calculated  according  to  equation  (1.2)  when  yq  =  1.55  g/cm^  and 
when  values  of  are  equal  to:  0.001;  0.01;  0.02;  0.03;  0.05;  0.1, 
respectively.  These  curves  determine  the  compressibility  of  the 
soil  as  a  three-component  medium. 

Fig.  4. 


The  stress  in  the  soil  in  the  direction  of  the  wave  movement 
is  equal  to  the  sum  of  P1  =  p-c^  stresses  in  the  skeleton  o.^  and 
in  the  soil  as  a  mixture  p.  In  the  case  of  unsaturated  sand  the 
curve,  calculated  according  to  equation  (1.2),  lies  much  lower 
than  curve  2,  p  <<  |o^|.  One  can  consider  therefore  that  curve  2, 
corresponding  to  the  overall  compressibility  of  the  soil,  practi¬ 
cally  coincides  with  the  curve  determined  by  the  compressibility 
of  the  skeleton.  The  compressibility  of  the  soil  as  a  mixture 
cannot  be  taken  into  account. 

With  the  change  in  the  content  of  water  in  the  pores  not 
only  does  the  compressibility  of  the  soil  change  as  a  mixture, 
but  the  compressibility  of  the  skeleton  of  the  soil  as  well,  since 
the  cohesive  forces  between  the  solid  particles  depend  on  the 
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moisture  content.  As  the  experiments  indicate  [3],  the  effect 
the  moisture  change  on  the  compressibility  of  the  skeleton  is 
substantially  less  than  the  compressibility  of  the  soil  as  a  mix¬ 
ture.  Therefore,  in  the  first  approximation  one  can  assume  that 
curve  2  determines  the  compressibility  of  the  skeleton  o1(e)  as 
well  as  the  water  saturated  soil.  Then,  curve  1  corresponds  to 
the  total  stress  in  water  saturated  soil  curve  2  -  to  the  stress 
in  the  skeleton  a^,  and  the  difference  of  these  stresses  —  the 
pressure  p  in  the  soil  as  a  mixture.  When  the  wave  is  not  a  shock 
wave,  curves  1  and  2  only  determine  the  general  character  of  the 
relationship  o1  =  o^Ce)  (dashed  lines). 

The  values  p  =  p1  +  are  close  to  the  stress  determined  by 
curve  6,  which  is  evidence  of  the  applicability  of  the  equation 
of  state  of  a  three-component  medium  (1.2)  to  the  soil  as  a  mix¬ 
ture.  Curve  6  corresponds  to  the  same  content  of  air  =  0.03 
as  there  would  be  in  a  saturated  soil. 

Let  us  examine  the  values  of  P1  and  p/P^  determined  by  curves 
1  and  2: 


P^^  kgf/cm2 

? 

p  kgf/cm 

P/P-L 

20 

15.6 

0.82 

14 

10. 

0.71 

8 

4.5 

0.56 

From  the  table  it  is  obvious  that  the  ratio  p/P-^  with  an 
increase  in  the  load,  increases  and  approaches  unity.  This  means 
that  with  the  larger  loads  the  compressibility  of  the  soil  to  large' 
degree  is  determined  by  its  compressibility  as  a  mixture  than  with 
smaller  loads. 

From  an  analysis  of  the  experimental  data  [3]  it  follows  that 
for  each  soil  the  value  of  the  load  =  p#  can  be  shown, 

beginning  whereby  one  can  ignore  the  compressibility  of  the 
skeleton  and  calculate  only  its  compressibility  as  a  three-component 
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medium.  With  an  increase  in  the  moisture  of  the  soil,  i.e.,  wit.. 

# 

a  lowering  of  a^,  the  values  of  are  reduced. 

* 

Approximate  values  of  P.^  for  water  saturated  soils  with  various 
cx^  are  given  below 

a1  <0.001  0.001-0.01  0.01-0.04 

P,  kgf/cm^  0  1-4  4-25 

* 

For  unsaturated  soils  the  values  of  P.,  are  equal  to  several 
hundreds  of  atmospheres. 

Within  the  examined  range  for  soils  of  high  and  low  moisture 
(Fig.  4)  we  will  have  correspondingly: 

P,  (*)  =  /»(*)-«,(«).  P,  (e)  =  -os(e).  (2.4) 

If  from  the  experiment  the  relationship  ?1  ~  P-^e)  is  known 
for  a  soil  of  high  moisture,  then  with  the  help  of  equations  (1.2) 
and  (2.4)  we  will  find  the  relationship  ox  ■  o1(e).  By  knowing 
from  the  experiment  the  relationship  P-^e)  *  -a^t)  for  a  soil  of 
low  moisture,  let  us  find  with  the  help  of  these  same  equations 
the  relationship  P1  =  P-L(e)  for  the  same  soil  at  high  moisture. 

3.  Let  us  examine  the  coefficient  of  lateral  pressure.  As 
the  experiments  indicate,  in  unsaturated  soil  when  one  can  take 
into  account  only  the  compressibility  of  the  skeleton,  ht  is  held 
practically  as  a  constant  value  with  an  increase  in  the  loads  up 
to  several  tens  of  atmospheres: 


k,  =  =  const  <  1.  (2.5) 


FTD-HT-2 3-7^9-72 
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In  a  water  saturated  soil  when  one  cannot  take  into  account 
the  compressibility  of  the  skeleton,  with  these  same  loads 

k*  - 1.  (2.6) 

In  a  soil  where  one  can  take  into  account  both  compressibility 
mechanisms,  the  stress  in  a  direction  perpendicular  to  the  direction 
of  the  external  load,  is  determined  by  the  equation 

Pa  (e)  =  p  (e)  -  Kx  o,  (•)•  (2.7) 

From  this  the  coefficient  of  lateral  pressure 


El.  =  P-wi 

Pt  p-a i 


(2.8) 


where  -  coefficient  of  lateral  pressure  determined  based  on  the 
calculation  of  the  compressibility  of  the  skeleton. 


If  the  external  load  is  not  too  small,  then  with  its  further 
increase  p  increases  to  a  greater  degree  than  does  a^.  From  (2.8) 
it  follows  that  in  this  case  the  ratio  P2/P-j_  increases  and  approaches 
unity.  The  properties  of  the  soil  approach  those  of  a  liquid. 


Let  us  present  the  experimental  values  of  P-^  and  ?2/'^i  a 
water  saturated  soil  which  confirms  this  deduction: 


P-^  kgf/cm£ 


P2/Pl 


7 

0.75 


15 

0.85 


20 

0.90 


The  average  normal 
pressibility  mechanisms 


stress  in  the  calculation  of  the  two  com 
is  determined  by  the  equation 
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3  =  -i-  (Pj  +  Pj  +  P|)  =  P - («I  +  8*  +  «»)• 


(-.9) 


In  this  way  a  model  of  the  soil  was  proposed  which  takes  into 
account  the  two  compressibility  mechanisms,  and  a  system  of  equations 
was  given  which  determines  the  planar  and  spherical  movements.  In 
the  case  of  non-unidimensional  movements  it  is  also  necessary  to 
take  into  account  the  two  compressibility  mechanisms.  In  this  case 
the  behavior  of  the  soil  as  a  mixture  is  determined  according  to 
the  equation  of  state  of  three-component  media,  proposed  by 
G.  M.  Lyakhova  [1],  and  the  behavior  of  the  same  soil  as  a  solid  - 
with  the  help  of  one  of  the  system  of  equations  proposed  in  [ h ,  5], 

The  authors  express  their  gratitude  to  R.  I.  Dubova  and 
V.  I.  Kuznetsov  for  their  participation  in  the  experiments. 
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